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We theoretically investigate the real-time emission spectrum of a two-level atom coupled to an 
optomechanical cavity (OMC). Using quantum trajectory approach we obtain the single-photon 
time-dependent spectrum in this hybrid system where the influence of a strong atom-cavity coupling 
and a strong optomechanical interaction are studied. We find a dressed state picture can explain 
the spectra by predicting the exact peak locations as well as the relative peak heights. In our 
analysis we also include the effect of mechanical losses (under weak mechanical damping limit) and 
single-photon loss through spontaneous emission from the two-level emitter. 


I. INTRODUCTION 


Hybrid atom-optomechanical quantum systems have 
shown a signihcant amount of research activity in last 
few years W- A mixture of two well developed helds 
of cavity quantum electrodynamics (CQED)[51 H] and 
quantum optomechanics [ZlIH], these hybrid systems are 
anticipated to be an integral part of future quantum tech¬ 
nologies where a wealth of phenomena can arise due to 
atom-light and light-mechanics interactions at the quan¬ 
tum level. Recently, it has been theoretically proposed 
and experimentally observed that these hybrid systems 
in which a mechanical oscillator interacts simultaneously 
with a radiation field and with a matter-like system 
(a single or a multi-atom system e.g. a Bose-Einstein 
condensate), can be used to implement many interesting 
effects of pure quantum nature. Some current and 
compelling examples include: Optomechanical cooling 
by A-configured atoms [9], entanglement generation and 
detection in atom-optomechnical systems m, quantum 
state control of Rydberg atoms using optomechanics m 
and generation of macro-level quantum superposition of 
nano-particles using spin-mechanical interactions in NV 
diamond centers inside nano-diamonds m- 
A relevant and important question in the context of 
such hybrid systems is what type of spectrum is emitted 
by such structures and what can be inferred from the 
spectrum about different features of atom-cavity and 
cavity-mechanics coherent interactions? In this paper 
we answer this question by theoretically exploring the 
single-photon time-dependent spectrum emitted by an 
atom-OMC system in the strong atom-cavity and strong 
optomechnaical interaction regimes (hereafter called 
the strong-strong coupling case). The experimental 
realization of a strong coupling at the single-photon 
level in a cavity QED setup has already been demon¬ 
strated in past m, while achieving single-photon strong 
optomechanical coupling in an OMC is anticipated to 
be realized in near future (as reflected by a substantial 
progress on the theoretical as well as experimental side 
in the single-photon optomechnaics [IdHlT] ). 

Recently we have investigated the time-dependent 
spectrum of a single-photon initially trapped in an 


empty OMC cavity [18]. Here we extend that study to a 
more interesting situation when the actual single-photon 
source (an excited two-level emitter) is coupled to the 
OMC. This problem (from the perspective of standard 
CQED setups) is more appealing as the observation 
of vacuum Rabi splitting [T5j in the single-photon 
CQED is a key indicator of achieving strong coupling 
between a two-level atom and an optical cavity [20] . 
In present study we explore the modifications strong 
optomechanical interaction brings into the Rabi-splitted 
spectrum and more importantly we would like to answer 
that how these modifications emerge in the spectrum 
as a function of time? Here it is worthwhile to mention 
that there are some studies conducted in recent past, in 
which a laser or single photon driven hybrid atom-OMC 
system’s spectrum is calculated in the stationary limit 
(time-independent version) pTII23| . But here we focus 
on the time-dependent/real-time version of the spec¬ 
trum, which to our knowledge has not be reported yet 
in the present context. 

We find that even working in the strong atom-cavity and 
strong optomechanical coupling regimes (which are two 
necessary conditions to observe vacuum Rabi splitting 
in the CQED setups and multiple side-band peaks in the 
context of cavity optomechanics, respectively), it takes 
finite time (t = 20u;]^) to observe the fully resolved 
spectrum in an atom-OMC system. Moreover, spectrum 
evolves from a broad Lorentzian to a frequency doublet 
and hnally additional side-band appears in the spectrum. 
Our analysis gives quantitative values of time when 
these spectral changes can be observed. The time order 
in which spectrum grows and manifests these features 
indicate different physical process responsible for the 
production of multiple resonances in the spectrum. 

The paper is organized as follows: In Sec. H we 
present our theoretical model describing the hybrid 
atom-optomechanical cavity setup. Next in Sec. HI 
we discuss our results in detail, focusing on the single¬ 
photon spectrum emitted in real time (time-dependent 
spectrum) when there is a strong coupling between atom 
and the cavity mode and also between single photon 
and movable cavity mirror in the OMC. A dressed state 
picture is employed to understand the spectrum by 
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exactly locating the spectrum resonances and relative 
peak heights. Furthermore the effect of phonon losses (in 
the waek mechanical damping limit) and the events of 
atomic spontaneous emission are also examined. Finally 
in Sec. IV we summarize the conclusions of this work. 


describes the atomic Hamiltonian, second term denote 
the atom-cavity coupling, third term the Hamiltonian 
of a free quantum mirror and finally the fifth term 
represents the non-linear optomechanical coupling. 


B. Quantum trajectory method 

II. THEORETICAL DESCRIPTION 


A. Model and Hamiltonian 

As shown in Fig. 1, system under consideration con¬ 
sists of an initially excited two level atom embedded in an 
optomechnical cavity. This two level atom has a ground 
state (I 5 )) and an excited state (|e)) with transition fre¬ 
quency and spontaneous emission rate given by Weg and 
7 a respectively. Experimentally such a two level system 
can be realized using a solid state quantum dot m or 
an NV center in diamond |24j or an artificial atom in the 
context of circuit quantum electrodynamics (CirQED) 
[25] . For the sake of simplicity atom’s external degrees 
of freedom are neglected and optomechnaical cavity is 
taken to be a usual Fabry-Perot cavity with a fixed but 
partially transmitting (movable but perfectly reflecting) 
mirror on left (right) side. Optical cavity is assumed to 
have a single isolated resonant mode given by frequency 
LOc while the annihilation of photons in the cavity is de¬ 
scribed by the action of operator a. The atom-field in¬ 
teraction strength is given by the rate ga- 
We have considered the movable mirror to be a quantum 
harmonic oscillator with equilibrium frequency lom and 
the destruction of quanta of vibration (phonons) is rep¬ 
resented in term of an operator b. The parameter gM is 
used to indicate the optomechanical interaction strength. 
The mechanical oscillator is coupled with a heat bath as 
well which causes phonons to decay at a rate 7 m and 
the bath temperature is described in terms of average 
phonon number M. The OMC can also leak out photons 
at a rate k into a dispersion-less optical fiber which is 
assumed to have a continuum of modes. Optical fiber 
guides the single-photon wavepacket such that finally it 
can be detected at the output detector D. Neglecting 
the optical and mechanical baths, the system Hamilto¬ 
nian under rotating wave approximation and in a frame 
rotating with frequency Wc can be expressed as: 

Hsys = -fiAaa-a+ + h{gaa^a- + g*aa+) + huMb^b 
— hgMa^aib^ -I- b) 

( 1 ) 

Here Aa = uieg — and we have neglected the zero 
point energies of the optical and mechanical oscilla¬ 
tors. Non-vanishing commutation relations are given 
by: [d-+,i 7 _] = (T 2 , [o, fll] = 1 and [b,b^ = 1. Here 
az = |e) (e| — \g) {g\ while a+la- are the atomic 
raising/lowering operators and Uz is the usual Pauli 
operator. On the right side of Eq. (1), the first term 


To study the dynamics of such a hybrid atom- 
optomechnaical system when it is coupled to fiber’s mode 
continuum, we employ the quantum trajectory/jump 
approach (QJA) [5SH2S]. This open quantum system 
method is known to be an appropriate technique as it 
incorporates both the presence of photo-detectors in the 
setup as well as different sources of photon and phonons 
losses. In QJA the registration of an excitation (either 
a photon or phonon) is described by a jump operator. 
In present situation we have three types of jumps pos¬ 
sible: optical jumps (Jout) which describe a situation in 
which a photon leaks out from the cavity into the fiber 
and detected at detector D, atomic jumps (J^at) through 
spontaneous emission mechanism and mechanical jumps 
Omt) representing phononic-decays. Theses stochasti¬ 
cally occurring jump events are described by the follow¬ 
ing set of operators: 


j(o) = v^a. 

( 2 a) 


( 2 b) 


( 2 c) 


According to QJA, in the absence of quantum jumps, 
system evolves according to a non-Unitary Schrodinger 
equation which is given by: 


dt 


Ht) 


Hnh 



(3) 


where Hnh is a non-Hermitian Hamiltonian constructed 
from the system Hamiltonian Hsys (which is Hermitian) 
and an anti-Hermitian Hamiltonian constructed from the 
jump operators as follows: 


Hnh = j(0)_^jt(A) j(A)_| jt(M) J(M) 

(4) 

Note that in writing above Hamiltonian we have assumed 
the mechanical heat bath to be at zero temperature (i.e. 


M = 0). The ket 


is the so called no-jump state 


which is constructed from all different possibilities of find¬ 
ing the excitations in the system before being leaked out. 
In present case no-jump state takes the following form: 


m(^) l^j 0; ^ ^ (^) 

m—0 m—0 


on the right hand side of above equation we are follow¬ 
ing the notational convention that the first slot in the 
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FIG. 1: (Color online) An initially excited two-level atom in an OMC. Initially the OMC is assumed to have 
zero photons and phonons. The optical field emitted by OMC is guided through an optical fiber to the detector 
D and its time-dependent spectrum is used to probe the signatures of atom-cavity and optomechanical 
interaction dynamics. For further deaths about the model of the system see Sec. 11(A). 


ket labels atomic state, second slot indicates the number 
of photons in the optical cavity and third slot represents 
the number of phonons in the mechanical oscillator. In 
present paper we’ll focus on a situation when there are 
neither photons nor phonons present in the system ini¬ 
tially i.e. initial no-jump state is |e, 0,0). Now by ap¬ 
plying the machinery of QJA, in the next section, we’ll 
calculate the time-dependent spectrum. 


C. Time-dependent spectrum 

To understand the growth of spectrum as a function of 
time, in this section we calculate the single-photon time- 
dependent spectrum produced by an atom-OMC system. 
For this calculation we utilize a setup in which the single 
photon wavepacket emitted by the atom-OMC systems 
passes through a filter cavity which has a fixed bandwidth 
F but a variable frequency u and then after leaving this 
filter cavity the photon is registered by detector D. The 
time-dependent spectrum is then interpreted as the fil¬ 
tered counting rate after many repetitions of the same 
experiment and it can be expressed by using the defini¬ 
tion first introduced by Eberly and Wodkiewicz (E&W) 
[29] as given by: 

^0 Jo 

{d\t )d{t ))dt dt . 

( 6 ) 

where A = uj — ujc and in the actual E&W spectrum def¬ 
inition we have replaced the classical field amplitudes by 
the corresponding optical field annihilation and creation 
operators. Sometime to make comparison with the al¬ 
ready reported stationary spectra, (as give by the Wiener 
Khinchin theorem [20|) above spectrum is also integrated 
over the entire time range. We’ll call this integrated 
spectrum as the synthesized spectrum Ns{t] A, F) and in 


present article we’ll always report this synthesized ver¬ 
sion of the spectrum. 


III. RESULTS AND DISCUSSION 

There are many parameters involved in the dynam¬ 
ics of the system under study, we here study the effect 
of varying three parameters on the single-photon time- 
dependent spectrum. These parameters represent dif¬ 
ferent loss mechanisms and are: cavity decay rate k, 
mechanical damping rate 7 m and spontaneous emission 
rate ja- We’ll start with the simplest situation when 
the atomic spontaneous emission and mechanical damp¬ 
ing are completely neglected (i.e. we are setting 7 m = 0 , 
M = 0 and 7 a = 0 for Sec. 111(A)). 


A. Absence of mechanical and spontaneous 
emission losses 


In this subsection we’ll assume that the cavity leakage 
rate is the main and dominant source of photon decay 
which can be realized experimentally by assuming that 
K ^ 7 a and phonons decay is completely ignored (i.e. 
7m = 0). Under these conditions time-dependent spec¬ 
trum takes the following form: 


7V(t;A,r) 



X bm{t')dt" 


2 

dt . 



^{iA+r/2+imuiM)t 


(7) 


In Fig. 2 we have plotted the single-photon time- 
dependent spectrum in the strong-strong coupling case 
(ga/n = 8, gn/n = 2.4). We have chosen this regime 
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FIG. 2: (Color online) Single-photon time-dependent spectrum in a hybrid atom-OMC system, calculated at six 
different times, in the strong-strong coupling regime and within good cavity limit, t = 2u!'^, Auj'^, 

7oj^, and 20u}'^ correspond to thick blue dashed, red dotted, green dotted dashed, thin magenta dashed, 

brown thin dotted dashed and black solid curves respectively. We have considered an on-resonance case i.e. 

Aa = 0 while other parameters used are: ga/t^M = 4, gn/^M = 1-2, = 0.5, T/ojm = 0.1. Notice that 

spectrum starts as a broad Lorentzian and then vacuum Rabi-spitting appears in the spectrum, finally blue and 
red side-bands emerge almost at the same time, which later on become more visible in the spectrum. In above 
plots, we have considered the possibility of having upto ten phonons in the mechanical motion. This choice is 
made noticing that with taking higher number of phonons (m > 10 ) the overall structure of the spectrum 
doesn’t change rather only more tiny side-bands appear at higher AA/lom values which lie outside the range of 
the A/wM-axis shown. 


because it is well known that under this choice of 
parameters an atom-cavity system and an empty OMC 
system, both separately show a fully resolved spectra 
[Min]- Note that in the case of OMC spectrum one 
has to respect the additional good cavity limit {lom > k) 
as well. 

Staring from Fig. 2(a), we notice that in the beginning 
(tcoM = 1) spectrum appears as a broad Lorentzian curve 
and it takes at least t = time for the spectrum 

to split into two peaks. This is an indication that the 
spectrum has started to exhibit vacuum Rabi splitting. 
But it takes further 2ui^j time for the spectrum to 
show full Rabi splitting i.e. both peaks in the strong 
atom-coupling regime get separated by a frequency 
equal to 2ga (as shown in Fig. 2(c) and as one expects 
from the standard Rabi CQED model). This specifies 
a physical process in which an initially excited atom 
emits a photon into the cavity mode which is exchanged 
between the atom and cavity few times before leaking 
out from the cavity (see Fig. 2(c)). Note that we have 
taken ga > givi thus even if the single photon undergoes 
Rabi oscillations before being leaked out from the cavity 


and produces full Rabi-splitting yet the optomechanical 
interaction doesn’t leave its strong signatures on the 
spectrum. 

Around t = we begin to notice optomechanical 

interaction imprints on the spectrum in the form of 
additional side-bands (tiny shoulders). These side-bands 
are showing that now the single-photon has interacted 
for long enough time with the movable mirror such that 
it can now produce observable effects in the spectrum. 
Note also that both red and blue side-bands are visible 
on the negative and positive A/wM-axis respectively. 
Physically red (blue) side-bands are produced when 
single photon interacts with the movable mirror and 
loses (gains) some energy from the mechanical oscillator, 
thus coming out from the cavity with lower (higher) 
frequency as discussed in detail in Ref. mm- 
In the case of single-photon time-dependent spectrum 
emitted form an empty OMC (as reported in our recent 
work m) the red side-bands should appear in the 
time-dependent spectrum before the emergence of blue 
side-bands. But here we notice that both side-bands 
appear simultaneously. This shows that in the process of 
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generating fully splitted vacuum Rabi spectrum photon 
has interacted with the mechanical mirror but due to 
small optomechanical interaction compared to Qa and 
due to broadening of vacuum Rabi spectrum peaks, red 
side-bands are not resolved before the blue ones. At 
later times (as shown in Fig. 2(e) at t = due to 

optomechnaical and atom-cavity interactions occurring 
for longer times, we hnd both of the Rabi peaks become 
sharper as well as both blue and red side-bands start 
to show up simultaneously at integer multiple of wm 
apart from the main Rabi peaks. Finally at t = 20uj^ 
we obtain the stationary limit of the spectrum as also 
reported in [22]. It is worthwhile to note that in Ref. 
[ 22 ] this stationary spectrum is calculated based on a 
real-space quantification approach as first introduced by 
Shen et.al in 2009 [3T]. Here, along with reporting the 
time-dependent version of the spectrum, we have also 
provided an alternate method of calculating the time- 
independent spectrum using the quantum trajectory 
method, which gives the exact same results as one can 
find by applying other methods for example the method 
of Shen et.al. 

In addition to the conclusions made in references [22] 
(to which our stationary spectrum results are in agree¬ 
ment), we conclude from the time-dependent spectrum 
that even if we are working in a strong-strong cou¬ 
pling regime still it takes ~ t = 2t)u}^ time to observe 
fully resolved spectrum. Moreover, the time-dependent 
spectrum also reveals the order in which different reso¬ 
nances appear in the spectrum describing the dynamic 
role played by different physical mechanisms which are 
responsible for the production of these spectral peaks. 
Finally at the end of this section, we would like to men¬ 
tion that by taking ga = 0 and qm = 0 in the stationary 
spectrum one can recover the usual single-photon spec¬ 


tra obtained in the context of CQED [2D] and an empty 
OMC [H] respectively. 


1. Dressed State Analysis and Asymmetry in peak heights 

In order to understand the stationary/time- 
independent spectrum we now present a dressed 
state analysis of the problem. For that we’ll consider a 
simple case when the atom-OMC starts off with atom 
in the excited state and movable mirror in its quantum 
ground state while later on there can be at most a single 
phonon present in the system. Note that this single 
phonon restriction is basically an artificial situation 
because the condition ^ << 1 should be satisfied for 

K 

the single-phonon approximation to be valid in an actual 
experiment. But the imposition of the condition — < 1 
is equivalent to working in the weak optomechanical 
coupling limit in which a fully resolved spectrum can’t 
be observed. Yet for simplicity (and to present analytic 
results) we’ll report this analysis here which can easily 
be extended (numerically) to as many as ten phonons 
case while in that situation the assumption of strong 
optomechnaical coupling is going to be a valid. 

To perform a dressed state analysis we neglect all 
source of decoherence and concentrate on the system 
Hamiltonian (as given in Eq. 1). With a single photon 
and a single phonon restriction one can define the ba¬ 
sis: {| 5 , 0 , 0 ), |g, 0 , 1 ), \g, 1 , 0 ) , \g, 1 , 1 ), |e, 0 , 0 ) , |e, 0 , 1 )} 
needed for the diagonalization of the system Hamilto¬ 
nian. As a result we obtain the following set of eigen¬ 
values (in unrotated frame) with different coefficients ap¬ 
pearing in the un-normalized eigenvectors shown in Fig. 3 
given by: 


Agoi,E±±| = jfiwg, fiwg-k fkUM, ^(wm ± ^gl + 2^^ -k ± 2 ^+ 9m + | 

-all + ± +5m + 45aWM ± ^M \J ± 2 ^Aglglj + gh + 4.glujlj 


a±±± = 


‘2gagM 


b±± = 


-WM ± y 4g2 -h 2glj +ujIj± 2^Aglgli + gij + 'igl 




‘^9a 


C(±±)(±±) 


-2(ga - (^^ ± + 2gM + cum ‘^V^^alah + 9 m + 


— 9a9M<^M 

2 


,2 


± 9agM\l‘^gl + ‘2gM + ± 2y/Aglglj + gh + 

wm ± ^(-WM ± ]J 4g2 -H 2glj + wir ± 2^Aglglj + glj + 


( 8 a) 

( 8 b) 

( 8 c) 

( 8 d) 


while in Fig. 3 we have E _= Ei,E _= E 2 ,E^ _= = E^. We note that due to atom-cavity and 






















6 



-leoi) 


-leoo) 


Bare/Uncoupled States 


Egii 

Egio 

Egoi 

Egoo 


■ Igll) 

■ IglO) 

Igoi) 

IgOO) 


[Dressed Stat^ 


E4 


E3 


E2 


(I) 


El . 
Egoi . 

EgOo ' 


(ii) 


a+-+lgio>+b-+lgii)+c(-+)(-+)leoo>+|eoi> 


(HI) 


(IV) 


ai-++lgio) + b++lgii>+c(++)(++)leoo>+|eoi> 


(V) 


a-lgio>+b-lgii>+C(++)(++)leoo>+|eoi> 

(VI) 

a lgio)+b+-lgii>+c(+-)(+-)l®°o>+ 


' i ' > 


(VII) 


(VIII) 

IgOl) 


IgOO) 


FIG. 3: (Color online) Dressed state picture of a hybrid atom-OMC system. Only single photon and single 
phonon restricted energy level diagram is shown. We notice that due to atom-cavity coupling and 
optomechanical interaction upper four states in the bare state picture are mixed up to form four dressed states. 
Transition from these excited states to two ground states are responsible for the production of eight resonances 
as shown in Fig. 4. Moreover this dressed state picture can also be used to explain the peak asymmetry present 
in Fig. 4 (see text for details). 


optomechanical interaction four bare states get mixed to 
form four dresses states. From the upper (excited) four 
dressed states system can make a photonic transition to 
two ground states. We note that this process gives rise to 
a total of eight transitions and these transition frequen¬ 
cies exactly matches the peak locations obtained in the 
spectrum as shown in Fig. 4. Here we would like to indi¬ 
cate that there exists an alternate method of calculating 
peak locations as well. This method is performed by set¬ 
ting the real part of the poles equal to zero in the final 
expression of the time-independent spectrum [ITj. We 
confirmed that the peak locations obtained from both 
methods (dressed state analysis and the pole-method) 
gives us identical results. 

It is worthwhile to mention here that the eigenvalues ob¬ 
tained upto single phonon are not in agreement with the 
result obtained in the reference [55] (energy eigenvalues 
vary as ~ bio'll -|- In this reference (and 

some other references mini related to the subject of 
single photon optomechanics) a polaron transformation is 
applied to transform the non-linear optomechanical cav¬ 
ity into a frame where the Hamiltonian can easily be dig- 
onalized. The reason for this disagreement again lies on 
our artificial single phonon restriction. We checked (nu¬ 
merically) that when more and more phonons are added 
in our model both results (one which is obtained by fol¬ 
lowing our quantum trajectory method and the other 
which performs a polaron transformation) start to match. 
Thus this issue can easily be resolved by extending the 
present single phonon restricted analysis to a multiple 



FIG. 4: (Color online) Single-photon single-phonon 
stationary spectrum emitted by an atom-OMC 
system. This spectrum is plotted to confirm that 
the dressed state analysis shown in Fig. 3 provides 
exact peak locations. Different peak labels 
corresponds to different transitions among dressed 
states as marked in Fig. 3. Parameters are the same 
as used in Fig. 2. 


phonon situation. 

Another advantage of performing a dressed state anal¬ 
ysis is that it also explains the asymmetry in the peak 
heights which is not possible to explain otherwise (for ex¬ 
ample by by setting the real part of pole equal to zero). 
In the dressed state picture we find that with the pres- 
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ence of single photon in the atom-OMC, system can be 
in any one of the four upper/excited states. In general, 
the state of the system at a later time t will be in a su¬ 
perposition of all four excited states and since the form 
of superposition can be imbalanced hence the transition 
from any one of these excited states to bottom ground 
states (I 5 OO), I 5 OI)) will happen with unequal probability 
thus causing an asymmetry in the peak heights. Finally 
we would like to remark that we numerically extended 
this dressed state analysis to the case of ten phonons 
case (corresponding spectrum shown in Fig. 2(f)) and we 
found that dressed state method gives exact peak loca¬ 
tions as well as describe the unequal peak heights for that 
case too (the analysis not shown here). 


IV. INCLUSION OF MECHANICAL DAMPING 
AND SPONTANEOUS EMISSION 


In this section we consider a more realistic situation 
in which two other sources of dissipation: mechanical 
damping and atomic spontaneous emission are also in¬ 
corporated in the model. To include mechanical decays 
we assume the movable mirror is coupled with a finite 
temperature Markovian mechanical heat bath. The tem¬ 
perature of the bath is characterized by the average ther¬ 
mal phonon number M and as a result of this coupling 
phonons can now leak out and enter into the OMC at 
a rate 7 m- Assuming the initial state of the mechanical 
oscillator to be in thermal equilibrium with the mechan¬ 
ical bath, we can express the initial thermal state of the 
oscillator using the Boltzmann probablity distribution as: 


Pm{t = 0) 


M ' 


(1 -bM)’"+i' 


(9) 


while label m is the number of phonons present in the 
movable mirror. Like previous section, we now apply 
the quantum trajectory method to calculate the time- 
dependent spectrum. Non-Hermitian Hamlitonian in this 
case takes the following form: 

^ JKO)J(O) - 
■n h 

_ |(M+ l)jt(M)J(M) _ 


with the presence of mechanical bath, the spectrum now 
become a weighted average over spectra calculated with 
different initial numbers of phonon states. Alongwith 
photonic detection events there will be undetected 
phononic quantum jumps occurring now as well which 
are represented by jump operator (see Eq. 2(c)). 


The jump and no-jump evolution leads to a net decay 
of phonons from the system which is given by the rate 
7m, as one can also see from the full Master equation 
analysis of the problem as well. In order to ensure that 
the optical measurement based state reduction to be a 


more probable situation than the mechanical decoher¬ 
ence in our setup, we’ll also assume a weak mechanical 
damping limit (M 1 such that k » 7 m-^) [5^ 155] . 
Here we would like to emphasize that although quantum 
trajectory theory can easily be applied to high phonon 
number case also, but since in present work we are more 
interested in observable photonic jumps (which will lead 
us to the optical spectrum) rather than the unobserved 
mechnaical jumps, hence we are considering a weak 
damping limit. Notice also that in high damping case, 
quantum trajectory theory will spend most of its time in 
calculating the unobserved mechanical quantum jumps, 
which we want to avoid in present optical spectrum 
study. This weak damping limit will allow us to neg lect 
the contribution coming from the M terms in Eq.(lO) 
compared to 1. 


The time-dependent spectrum incorporating phononic 
and atomic decays is plotted in Fig. 5. Like Fig. 2 we ob¬ 
serve that in order to resolve the spectrum into vacumm 
Rabi spitting it takes some finite time (~ t = as 

shown in Fig. 5(c)). The novel part of the spectrum is 
that now the phononic side-bands takes even longer to 
show up in the spectrum. In the absence of 70 and 7 m 
we found that the mechanical motion signatures start to 
appear in the spectrum much earlier (~ t = as 

shown in Fig. 2(d)), but when phonon losses are taken 
into account we observe the emergence of phononic side¬ 
bands happens later around t = IQu)^. 

This behavior can be explained based on the fact that it 
takes some finite time for the single-photon to interact 
with the movable mirror and to carry the information 
about this interaction in its spectrum. When there is no 
mechanical damping present phonons generated due to 
optomechanical coupling will stay in the OMC and hence 
photon will exhibit the phonon presence fastly. Com¬ 
pared to that, when phonons can be lost from the OMC 
quickly, then photon will need more time to show ob¬ 
servable effects produced by photon-phonon interaction 
(phononic side-bands), as one can notice by comparing 
Fig. 2(d)-(e) and Fig. 5(d)-(e). Similar to Fig. 2, both 
red and blue side-bands appear at the same time in the 
spectrum now. 

Finally at t = (Fig- 2 (f)) we obtain the infinite 

long time version of the spectrum. We note that al¬ 
though the spectrum has reduced its height considerably, 
still major vacuum Rabi splitted peaks as well as signs of 
side-bands are visible. We found also that both mechan¬ 
ical damping and spontaneous emission contribute in the 
overall peak reduction compared to Fig. 2(f). Through 
spontaneous emission process this happens due a de¬ 
crease in the overall photon detection events as consid¬ 
erable number of times the emitted photon is lost into 
unwanted spontaneous emission channel (a process hap¬ 
pening at a rate 7 a) and hence filtered counted rate (time- 
dependent spectrum) as displayed in Fig. 5(f) shows con¬ 
siderable height reduction compared to Fig. 2(f). Me¬ 
chanical damping not only causes peak height reduction 
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FIG. 5: (Color online) Single-photon time-dependent spectrum in a hybrid atom-OMC system, including 
mechanical damping and spontaneous emission processes. The spectrum is recorded as six different times. 
t = ,2uj]^ and 20uj'^ which correspond to thick blue dashed, red dotted, green dotted 

dashed, thin magenta dashed, brown thin dotted dashed and black solid curves respectively. The parameters 
used are: 7 m/wm = 0.1, ya/wM = 0.4, M = 0.1 and rest of the parameters are the same as used in Fig. 2. Note 
that we have chosen weak damping but a large spontaneous emission loss {ja k) in plotting this figure. 


but also causes a broadening of peaks due to blurring of 
resonances in the wide thermal background. Note that 
such type of mechanical blurring effect is also discussed 
in detail in Ref. [16], where an empty OMC is driven 
by a laser field and a cross over from multi-photon to 
single-photon optomechanics is investigated. 

V. CONCLUSIONS 

By employing the quantum jump theory, we have stud¬ 
ied the single-photon time-dependent spectrum emitted 
from a hybrid atom-OMC system in a strong-strong cou¬ 
pling regime. We concluded that the time-dependent 
spectrum shows how the spectrum builds up in time 


staring form a broad Lorentzian to a vacuum Rabi split- 
ted frequency doublet and finally showing the emergence 
of phononic side-bands. Our quantitative analysis pro¬ 
vides exact times when optical spectrum suffers these 
changes. The order in which these resonances appear 
also describe the time order in which different physical 
processes become dominant in the system. We found a 
dressed state picture can be utilized not only to predict 
the exact locations of different peaks appearing in the 
spectrum but also to understand unequal peak heights. 
Finally, by including mechanical losses and spontaneous 
emission events, we found that heights of all peaks ap¬ 
pearing in the spectrum show marked decrease such that 
tiny side-bands appearing at large zLA/ujm values dimin¬ 
ish completely. 
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